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Abstract 

The Jones polynomial, discovered in 1984 I18| . is an important knot invariant in topology. Among 
its many connections to various mathematical and physical areas, it is known (due to Witten P2j ) to be 
intimately connected to Topological Quantum Field Theory (TQFT). The works of Freedman, Kitaev, 
Larsen and Wang I13II14I provide an efficient simulation of TQFT by a quantum computer, and vice versa. 
These results implicitly imply the existence of an efficient quantum algorithm that provides a certain 
additive approximation of the Jones polynomial at the fifth root of unity, e^'^'''^, and moreover, that 
this problem is BQP-complete. Unfortunately, this important algorithm was never explicitly formulated. 
Moreover, the results in |13l I14| are heavily based on TQFT, which makes the algorithm essentially 
inaccessible to computer scientists. 

We provide an explicit and simple polynomial quantum algorithm to approximate the Jones poly- 
nomial of an n strands braid with m crossings at any primitive root of unity g^'^V*:^ where the running 
time of the algorithm is polynomial in m,n and k. Our algorithm is based, rather than on TQFT, on 
well known mathematical results (specifically, the path model representation of the braid group and the 
uniqueness of the Markov trace for the Temperly Lieb algebra). By the results of |14|. our algorithm 
solves a BQP complete problem. 

The algorithm we provide exhibits a structure which we hope is generalizable to other quantum 
algorithmic problems. Candidates of particular interest are the approximations of other downwards self- 
reducible T^P-hard problems, most notably, the important open problem of efficient approximation of 
the partition function of the Potts model, a model which is known to be tightly connected to the Jones 
polynomial |33| . 

1 Introduction 



Since Shor's breakthrough discovery in 1994 |28| . quantum algorithms with an exponential speedup over 
the best known classical algorithms have been shown for a number of problems (e.g., |10[ I30L [TBI \2'2\ ). All 
these problems and algorithms share some common features: the problems are group or number theoretic in 
nature and the key component of each algorithm is the quantum Fourier transform^. Arguably, the greatest 
challenge of quantum computation is the discovery of new algorithmic techniques. 

In this paper we describe a polynomial time quantum algorithm that approximates the ^P-hard problem 
of evaluating the Jones polynomial at certain roots of unity. The best classical algorithm for this problem 
is exponential. Our algorithm is significantly different from all previously known quantum algorithms that 
achieve an exponential speed up in the following three important ways: 1) it solves a problem which is 
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combinatorial rather than group or number theoretic in nature, 2) it does so not by using the Fourier 
transform, but instead, by exploiting a certain structure of the problem and encoding it into the nature 
of the unitary gates being used, and 3) it solves a problem that is BQP-hard ^21j that is, a problem that 
captures all the power of the quantum model. 

The connection between quantum computation and the Jones polynomial was first made through the 
series of papers jlll 1121 1131 114| . A model of quantum computation based on Topological Quantum Field 
Theory (TQFT) and Chern-Simons theory was defined in |lllll2j . and Kitaev, Larsen, Freedman and Wang 
showed that this model is polynomially equivalent in computational power to the standard quantum com- 
putation model in '13[ I14|. These results, combined with a deep connection between TQFT and the value 
of the Jones polynomial at particular roots of unity discovered by Witten 20 years ago |32| , implicitly 
implies an efficient quantum algorithm for the approximation of the Jones polynomial at the value e^^'/^. 
This connection is also discussed, from the point of view of TQFT, in Preskill's notes |2S]- Unfortunately, 
the important quantum algorithm implied by these intriguing results, though referred to in \7\, was never 
explicitly formulated. 

In this paper we use a different route to connect quantum computation and the Jones polynomial, one 
that does not involve TQFT. We present an explicit and simple to state algorithm for the above problem, 
which is based purely on algebraic results from more than 20 years ago. Moreover, our algorithm works for 
all roots of unity of the form e^'^*/*'', going beyond the discussions in previous works involving only constant 
k's. We now describe the precise problem that we solve. 

1.1 Background on the Jones Polynomial 

A central issue in low dimensional topology is that of knot invariants. A knot invariant is a function on knots 
(or links -i.e. circles embedded in R^) which is invariant under isotopy of the link, i.e., it does not change 
under stretching, moving, etc., but no cutting. In 1984, Jones 18- discovered a new knot invariant, now 
called the Jones polynomial VL{t), which is a Laurent polynomial in Vt with integer coefficients, and which 
is an invariant of the link L. In addition to the important role it has played in low dimensional topology, 
the Jones polynomial has found applications in numerous fields, from DNA recombination |24j . to statistical 
physics Plj. 

The Jones polynomial can also be defined as a function of braids. A braid of n strands and m crossings 
is described pictorially by n strands hanging alongside each other, with m crossings, each of two adjacent 
strands. A braid may be "closed" to form a link by tying its ends together. In this paper we will be interested 
in two ways to perform such closures, namely, the trace closure and the plat closure (to be defined later). 
We will be interested in the Jones polynomial of links that are trace or plat closures of braids. 

The essential aspect of the Jones polynomial of a link L can be computed by the following process: project 

L to the plane keeping track crossings to get what is called a link diagram. Now replace every crossing \ 



by a particular linear combination (with coefficients being functions of the parameter t) of the pictures 



diagrams with a particular function of the parameter t and the number of loops. The resulting function of t 
is a scaled version of the Jones polynomial Vl (t) . 

From the moment of the discovery of the Jones polynomial, the question of how hard it is to compute 
was important. There is a very simple inductive algorithm (essentially due to Conway [3]) to compute it 
by changing crossings in a link diagram, but, naively applied, this takes exponential time in the number of 
crossings. It was shown that the computation of Viit) is ^j^P-hard for all but a few values of t where 
VL{t) has an elementary interpretation. Thus a polynomial time algorithm for computing VlII) for any value 
of t other than those elementary ones is unlikely. Of course, the ^P-hardness of the problem does not rule 




and 



. The result is a linear combination of diagrams containing only closed loops. Replace each of these 
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out the possibility of good approximations; see, e.g., ^Jj. Still, the best classical algorithms to approximate 
the Jones polynomial at all but trivial values are exponential. 



1.2 Our Results 

We show an efficient, explicit, and simple quantum algorithm to approximate the Jones polynomial at all 
points of the form t — g^in/fc, jt will in fact be easier to use the parameter A = t^^^'^. We prove the following 
for the trace closure B*'' and the plat closure -B^' of a braid B: 

Theorem 1.1 For a given braid B with n strands and m crossings, and a given integer k, there is a quantum 
algorithm which is polynomial in n, m, k which with all hut exponentially small probability, outputs a complex 
number r with \t — Vgtr (e^'^'/'^)| < ed^^^ where d = —A^ — A^"^ , and e is inverse polynomial in n, k,m. 

Theorem 1.2 For a given braid B with n strands and m crossings, and a given integer k, there is a quantum 
algorithm which is polynomial in n, m, k which with all but exponentially small probability, outputs a complex 
number r with \r — Vbpi (e^'^'/'^)| < ed^'^/'^/N where d = —A^ — A^^ and e is inverse polynomial in n,k,m 
(N is an exponentially large factor to be defined later). 

We remark that the approximation we provide here is additive, namely the result lies in a given window, 
whose size is independent of the actual value we are trying to approximate. This of course is much weaker 
than a multiplicative approximation, which is what one might desire (see discussion in One might 

wonder if under such weak requirements, the problem remains meaningful at all. It turns out that, in fact, 
this additive approximation problem is hard for quantum computation: 

Theorem 1.3 Adapted from Preedman, Larsen and Wang |14j The problem of approximating the 
Jones polynomial of the plat closure of a braid at e^'^*/'^' for constant k, to within the accuracy given in 
Theorem \l.Sl is BQP -hard. 

This result was recently strengthened by Aharonov and Arad ^ to any k which is polynomial in the size 
of the input, namely, for all the plat closure cases for which our algorithm is polynomial in the size of the 
braid. 

Curiously, the hardness results of ^^Q] are not known to hold (regardless of fc) for the approximation of 
the trace closure for which we give an algorithm as well. We discuss the difference between the two problems 
in the open questions section. 



1.3 Description of the Algorithm 

The essence of the algorithm lies in the fact that for braids with n strands, the pictures that one gets by 

replacing each crossing by one of the two pictures { ^ , }, can be assigned the structure of an algebra. 

This algebra is called the Temperley Lieb algebra, and denoted by TL^. In fact, the map from the crossing 
to the appropriate linear combination of the above two pictures defines a representation of the group i?„ of 
braids of n strands, inside the TL„ algebra. The Jones polynomial of the trace closure of a braid can be seen 
as a certain trace function (i.e., a linear function satisfying tr{AB) — tr(BA)) on the image of the braid in 
the TLn algebra. 

Our goal is then to design an algorithm that approximates this trace. To this end we use an important 
fact about this trace: it satisfies an additional property called the Markov property. Moreover, this property 
makes it unique; any trace function on the TL„ algebra (or a representation of it) that satisfies this property 
is equal to the above trace! This leads us to the key idea of the algorithm: suppose we can define a 
representation of the TL„ algebra by matrices operating on qubits, and we can identify and estimate the 
trace that satisfies the Markov property on this representation. Then by the uniqueness of this trace we can 
estimate the Jones polynomial. 
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But what is the representation that should be used? If our intent is to design a quantum algorithm, it 
is best if the representation induced on the braid group be unitary, so that we can hope to approximate its 
trace by a quantum computer. Fortunately, it is in fact possible to give representations of the Temperley 
Lieb algebra which induce unitary representations of the braid group. Such representations were constructed 
in |20[I18| and are called the path model representations. If we want to evaluate the Jones polynomial VL{t) 
for L a closure of a braid in i?„, and t — e^'^'/'^, we use the fcth path model representation of i3„. It is fairly 
straight forward to adapt these representations to work on the space of n qubits, and moreover, to show 
that the resulting unitary matrices (namely, the images of the generators of the braid group) can be applied 
efficiently by a quantum computer. We find that the image (by the path model representation) of an entire 
braid B can be applied efficiently by a quantum computer. Let us call the unitary matrix corresponding to 
a braid Q{B). 

To approximate the Jones polynomial of a trace closure of the braid, it suffices to approximate the Markov 
trace of Q{B). This is done using standard quantum and classical algorithmic techniques, including the well 
known Hadamard test. The algorithm for the plat closure builds on similar ideas, though it is not directly 
stated in terms of traces. Thus we obtain a polynomial quantum algorithm for the BQP-complete problem 
of approximating the Jones polynomial of a plat closure of a braid. 

We remark that after we completed this work, we learned about a previous independent attempt to prove 
similar results [201 • Unfortunately, the work of [221 is greatly flawed, and in particular claims to provide an 
exact solution to the ^P-\iaid problem. 

1.4 Conclusions and Further Directions 

We have provided a simple algorithm for a BQP-complete problem, which is different in its methods than 
previous quantum algorithms. In essence, what it does is to isolate a certain local structure of the problem, 
and assign gates which somehow exhibit the same local structure. Our hope is that this more combinatorial 
direction in quantum algorithms will lead to further progress in the area. 

In particular, one very interesting related question is an important problem from mathematical physics: 
that of approximating the partition function of the Potts model \M\ , which is known to be tightly connected 
to the Jones polynomial i and its exact evaluation is once again #P-hard [2] ■ We hope that the results of 
this paper will lead to progress in this question, or in other questions related to approximating #P-complete 
problems. 

We briefly discuss the relation between the plat and the trace closures problems. It is known that any 
plat closure of a braid can be transformed efficiently into a trace closure of some other braid 29 . The reader 
might therefore find it curious that one of these problems is BQP-complete, while the other one is not known 
to be so. The explanation lies in the fact that the quality of the approximation in both algorithms depends 
exponentially on the number of strands in the braid. The transformation from plat to trace closures requires, 
in the worst case, a significant increase in the number of strands. This, unfortunately, degrades the quality 
of the approximation exponentially. The computational complexity of the trace closure problem is left open. 

Finally, we believe that this paper helps to clarify and demystify (at least one direction of) the intriguing 
equivalence between quantum computation and the problem of approximating the Jones polynomial. We 
hope this connection leads to a deeper understanding of quantum computation complexity. 

Organization of paper: Section [3 provides the background and the necessary definitions, starting from 
quantum computation, the Hadamard test, the braid group, algebras and representations, the Jones polyno- 
mial, the Temperley Lieb algebras and the path model representation. Using these notions we describe the 
algorithms in Section |31 and prove their correctness. 
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2 Background 



2.1 Quantum Computation 

For background on quantum computation, see 123! . We merely mention here that, strictly speaking, we use 
here a quantum-classical hybrid model of computation, in which a classical probabilistic Turing machine 
performs calls to a quantum computer, and uses its outcomes to perform some classical computations. It is 
standard that this model can be simulated efficiently by the standard quantum computation model. 



The following fact is standard in quantum computation. If a state \a) can be generated efficiently, and a 
unitary Q can be applied efhciently, then there exists an efficient quantum circuit whose output is a random 
variable e {—1, 1}, and whose expectation is TZe{a\Q\a) . Start with the two-register state :^(|0) -I- 
apply Q conditioned on the first qubit to get the state :^(|0) |a) -I- |1) ® Q\a), apply a Hadamard gate 
on the first qubit, and measure. Output 1 if the measurement result is |0), —1 if the measurement result is 
1 1). The expectation of the output is exactly TZe{a\Q\a) . To get a random variable whose expectation is the 
imaginary part, start with the state -^(|0) — ® \a) instead. 

2.3 Algebra Background 

An algebra is a vector space with a multiplication. The multiplication must be associative and distributive. 
A representation of a group G inside an algebra is a group homomorphism p from G to the group of invertible 
elements in the algebra, namely, we require p{gi)p{g2) — p(.9i52) for any (71,(72 G G. 

We say a representation is reducible if there exists a proper subspace of vectors which is invariant under 
the group action. If there is no such subspace, we say the representation is irreducible. 

We shall sometimes refer to a representation of G without specifying the algebra; in these cases we shall 
mean a representation inside the algebra of n x n matrices. We shall be interested in algebra representations 
as well: 

Definition 2.1 An r dimensional representation $ of an algebra is a linear mapping from the algebra into 
the set ofrxr complex matrices Mr, such that for any two elements X, Y in the algebra, ^{X)^{Y) = ^{XY) 

If a group is represented inside an algebra then any representation of the algebra gives a representation 
of the group by composition. 

Often, an algebra or a group is defined using a set of generators and relations between them. In this case, 
a representation may be defined by specifying the images of the generators, provided the same relations hold 
for the images as for the generators. 



Consider two horizontal bars, one on top of the other, with n pegs on each. By an n strand braid we shall 
mean a set of n strands such that: (1) Each strand is tied to one peg on the top bar and one peg on the 
bottom bar, (2) Every peg has exactly one end attached to it, (3) The strands may pass over and under each 
other, (4) The tangent vector of every strand at any point along the path from top to bottom always has a 
non-zero component in the downward direction. Here is an example of a 4-strand braid: 



2.2 The Hadamard Test 



2.4 The Braid Group 
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The set of n-strand braids, Bn, has a group structure with multiphcation as follows. Given two n-strand 
braids 61, 62, place braid bi above &2, remove the bottom 61 bar and the top 62 bar and fuse the bottom of 
the 61 strands to the top of the 62 strands. 



The product of the above 4-strand braid with the 4-strand braid 



is: 





An algebraic presentation of the braid group due to Artin is as follows [3]: Let i3„ be the group with 
generators {1, cti, . . . (t„^i} and relations 

1. (TiCTj ~ (TjCJi for \i — i\ > 2, 

2. <7iat+i(7i ~ (Ti+icr.j(Ti+i. 

This algebraic description corresponds to the pictorial picture of braids: Oi corresponds to the pictorial 



braid 



and concatenating such pictures gives a general braid in i?„. 



2.5 The Temperley-Lieb Algebras 

Definition 2.2 Given n an integer and d a complex number we define the Temperley-Lieb algebra TL„{d) 
to be the algebra generated by {1, Ei, . . . , En-i} with relations 

1. E,E,^E,E,, \i-j\ > 2, 

2. EiEi±iEi = Ei, 

3. Ef = dE^. 

There is a well known geometric description of TLn{d) due to Kauffman It uses the notion of 

Kauffman n-diagrams, which is best explained by an example, e.g., a Kauffman 4-diagram: 




In general, a Kauffman n-diagram is a diagram as above, with n top pegs and n bottom pegs, and no 
crossings and no loops. More formally: 

Definition 2.3 Let Dn be a rectangle with n marked points on the top of the boundary and n marked points 
on the bottom. A Kauffman n-diagram is a picture sitting inside Dn consisting of n non-intersecting curves 
that begin and end at distinct marked boundary points. We will consider two such diagrams equal if they are 
isotopically equivalent (keeping the boundary fixed). 
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We define a vector space over these diagrams: 



Definition 2.4 Let ICn be the vector space of complex linear combinations of Kauffman n-diagrams. 

Multiplication of two Kauffman n-diagrams is done just like in the case of braids. To multiply a diagram 
ki with a diagram k2 we stack ki on top of fc2, and fuse the matching ends of the strands. The resulting 
diagram is a new Kauffman diagram fca with possibly some extra closed loops. We define the product of ki 
and k2 to be (F^kj G JCn where m is the number of extra closed loops. For example, a multiplication of the 

above Kauffman 4-diagram with the Kauffman 4-diagram results in: 





This multiplication rule can be extended linearly to /C„, we call the resulting algebra gTLn{d). 
The algebras TLn{d) and gTLn{d) are isomorphic: 



Theorem 2.1 The map ip : TLn{d) — > gTLn{d) given by the homomorphic extension of'ip{Ei) = 



is an isomorphism. 



17" 



Proof: It is a simple and fun exercise to check that the image of the relations given in Definition 12.21 are 
relations in gTLn{d). For the remaining details see □ 

We shall refer to the pictures of the form ^ , which generate gTLn{d), as capcups. 



2.6 Representing _B„ Inside TLn{d) 

We define a mapping from the braid group to TLn{d): 

Definition 2.5 pA '■ Bn TLn(d) is defined by 

PA{a,) = AE, + A-^I. 

Claim 2.1 For a complex number A which satisfies d = —A^ — A^"^ , the mapping pA is a representation of 
the braid group i3„ inside TL„{d). 

Proof: We need to check that the relations of the braid group are satisfied by this mapping. For \i — 
j\ > 1, pAio'i) commutes with pA{crj) since Ei commutes with Ej. To show that pAici) pAiui+i) pAio'i) = 
PA{o'i+i)pA{o'i)pA{o'i+i), substitute to get an expression in Ei's. Opening up the first expression we get 
A^EiEi+iE,, + AEi+iEi + AEf + A^^Ei + AE^Ei+i + A^^Ei+i + A^'^Ei + A^^. The second expression gives 
A^Ei+iEiE^+i+AEiEi+i+AEf_^_j^+A^^Ei+i+AEi+iE^+A~^Ei+A^^Ei+i+A^^. We remove similar terms, 
and using the relations of the TLn{d) it remains to show that {A"^ + Ad + A^)Ei — {A^^ + Ad + A^)Ei+i. 
This holds because the constants are due to the relation between d and A. □ 
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2.7 Unitary Representation of Bn 

Given a representation r of TLn{d), we may use the representation of the braid group inside the TLn{d) 
algebra (Definition 12. 5|l to derive a representation of Bn by composition, as follows. Define the map Lp by 
specifying its operation on the generators at of i?„ to be f{cfi) = tfi = T{pA{cri)) — AT{Ei) + A~^I. This 
representation is unitary under certain constraints: 

Claim 2.2 // \ A\ = 1 and T{Ej) are Hermitian for all i, then the map ip is a unitary representation of B^. 

Proof: T{pAia,))TipA{a^)V = {A'' I+At{E,)){{A-^)* I+A*TiE,y) = I+A-ME^)+AME^)+dT{E,) - /. 
□ 

2.8 Tangles 

For this paper, we define a tangle to be a braid in which some of its crossings have been replaced by a picture 
of the form { ^ }. Braids and Kauffman diagrams are tangles. 

2.9 From Braids to Links 

We can connect up the endpoints of a braid in a variety of ways to get links. We single out two such ways: 

Definition 2.6 The trace closure of a braid B shall be the link achieved by connecting the strand at the 
rightmost top peg, around to the right, to the strand at the rightmost bottom peg, then connecting in the same 
way the next to rightmost top and bottom strands and so on. We denote the resulting link by B*'' . 

Definition 2.7 The plat closure of a 2n-strand braid shall be the link formed by connecting pairs of adjacent 
strands (beginning at the leftmost strand), on both the top and bottom. We denote the resulting link by B^^ . 

Examples of the trace closure and the plat closure of the same 4-strand braid are: 



These closures are also well defined for tangles. 
2.10 The Jones Polynomial 

A definition of the Jones polynomial VL{t) due to Kauffman [2] is as follows. We start by defining the 
Kauffman bracket < L > as a polynomial in A for A such that = t. 

Definition 2.8 Consider a link L, given by a link diagram. A state a of L shall mean a choice, at each 
crossing ^ of L, from the set { ^ }. To a state a of a link L we associate the following expression 



<y{L): Let (respectively a ) he the number of crossings for which a chooses ^ (respectively 



)■ 



Let \a\ be the number of closed loops in the diagram gotten by replacing each crossing \ by the choice 
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indicated by the state a. Define a{V) = A'^"^ ^ . The Kauffman bracket polynomial, also called the 

bracket state sum, for L, is defined to be 

<L>= Yl ^(^)- 

all states a 

To define the Jones polynomial, we consider oriented links, namely links with one arrow on each connected 
component. The connection between the Jones polynomial and the Kauffman bracket is given by a notion 
called the writhe: 

\/ 

Definition 2.9 For an oriented link L, assign to each crossing that looks like this the value +1, and 

\/ 

to each crossing that looks like this: the value —1. The writhe of L is the sum over all the crossings 

of these signs. 

Definition 2.10 The Jones polynomial of an oriented link L is defined to be 

VL{t) = Vl{A-^) = <L> 

where w{L) is the writhe of the oriented link L, and < L > is the bracket state sum of the link L, ignoring 
the orientation. 

Thus, the Jones polynomial is a scaled version of the bracket polynomial. Moreover, the writhe of a link can 

be easily calculated from the link diagram, and hence the problem of calculating the bracket sum polynomial 
is equivalent in complexity to that of calculating the Jones polynomial. 



2.11 The Markov Trace 

Definition 2.11 A linear function from an algebra to the complex numbers is called a trace if it satisfies 
tr{XY) = tr{XY) for every two elements X, Y in the algebra. 

We define the following trace on gTL^id). 

Definition 2.12 The Markov trace tr : gTLn{d) C is defined on a Kauffman n-diagram K as follows. 

Connect the top n labeled points to the bottom n labeled points of K with non-intersecting curves, as in the 
trace closure. Let a be the number of loops of the resulting diagram. Define tr{K) = Extend tr to all 

of gTLnid) by linearity. 

For example: 




Since TLn{d) and gTL^id) are isomorphic, tr induces a trace on TLn{d); for simplicity we shall denote 
this map by tr as well. 
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Claim 2.3 tr satisfies the following three properties: 

1. tr{l) = 1, 

2. tr{XY) = tr{YX) for any X, F G TLr,{d), 

3. IfX e TLn-iid) then tr{XEn-i) = jtriX). 

Proof: It is straightforward to verify this by examining the appropriate pictures in gTLn{d). □ 

Of particular importance is the third property, which is referred to as the Markov property. These three 
properties uniquely determine a linear map on TLn{d): 

Lemma 2.1 Mfi^ There is a unique linear function tr on TLn{d) (and on any representation of it) that 
satisfies properties 1 — 3. 

Proof: By a reduced word w G TLn{d) we shall mean a word in the set {1, -Bi, . . . En-i} that is not equal 
to cw' for any c a constant and w' a word of smaller length. Using the relations of TLn{d) and applying 
simple combinatorial arguments, we show that a reduced word w G TLn{d) contains at most one En-i term. 
We induct on n. Clearly the only reduced words in Tl2{d) are 1 and Ei. Assume the statement is true for 
reduced words in Tln-i{d). Suppose there exists a reduced word w G TLn{d) containing more than one 
En~i term. Write w = wiEn-iW2En-iW3 with W2 a word without Since W2 must be reduced and is 

in Tl„-i{d), the induction hypothesis implies W2 contains at most one En-2 term. If W2 does not contain a 
En-2 term, W2 G Tln^2{d) and it commutes with so we have w — wiW2En-iEn^iW3 which shows that 

w was not reduced. Otherwise we can write W2 — vEn-2v' with v,v' both words in TLn-2{d). It follows 
therefore that v and v' commute with and thus w — wiv En-iEn-2En-iv' W3 which again shows that 

w was not reduced. We conclude that any reduced word in TLn{d) contains at most one En-i term. 

Given w G TLn{d)\Tln-i{d) a reduced word we write w — wiEn-iW2 with wi, W2 G Tln-i{d). Then 
tr(w) = tr{w2WiEn-i) — dtr{w2Wi), the first equality by property 2, the second by property 3. Thus for any 
word w G TLn{d) we can reduce the trace computation to the trace of a word W2W1 G TZ„_i(d). Iterating 
this process, (and using the fact that tr{l) = 1), we see that the trace of a word in TLn{d) is uniquely 
determined by the relations 1.-3. Since the trace is linear, the result follows. □ 

We have the following convenient description of the Jones polynomial in terms of the Markov trace. 

Lemma 2.2 Given a braid B, then 

Vb^^{A-^) = (-A)3»(^"-)d"-4r(p^(i?)). 

Proof: By Definition 12. 101 we need to show that < B*^ >= tr{pA{B))d^^^ . We observe that there exists 
a one to one correspondence between states that appear in the bracket sum < B*^ >, and Kauffman n- 
diagrams that appear in pa{B). The weight of an element in the bracket state sum corresponding to the 
state a is A°' d'""'^^. We observe that the corresponding Kauffman n-diagram appears in pa{B) with 
the weight A'^ ~°' . Hence, by linearity of the trace, it remains to show that for each a, the trace of the 
Kauffman diagram corresponding to cr, times d"^^, equals to the remaining factor in the contribution of a 
to the bracket state sum, d''"!"^. This is true since by the definition of the trace of a Kauffman diagram, it 
is exactly dl'^l"". □ 

This lemma also holds if B is replaced by a tangle. 

2.12 The Path Model Representation of TL^id) 

We describe the path model representation of TLn (d) due to |20l [H^ . The representation will act on a 
vector space determined by paths on a graph. Specifically, given an integer k [k will be chosen in relation 
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to d later), let Gk be the straight line graph with k — 2 segments and k — 1 vertices: 



« * • « • 




k-1 vertices 



Define Qn,k to be the set of all paths of length n on the graph Gk beginning at the leftmost vertex. Given 
q € Qn,k, we shall denote by ^(0), g(l), . . . q{n) the sequence of vertices of Gk describing q; thus q(0) is the 
leftmost vertex and q{i) and q{i + 1) are adjacent vertices of Gk for all i. We shall think of the elements of 
Qn,k as an orthonormal basis of a vector space Vn.k'i hence an element q € Qn,k shall represent both a path on 
Gk and a basis element of Vn,k- We now construct the path model representation T{TLn{d)) : Vn,k ^ ^n,k- 

Given a Kauffman n-diagram T, to describe r(T) it will suffice to give the matrix entry T(T)g/ g for each 
pair q',q & Qn,k- To do this, we note that the strands of a Kauffman diagram separate the rectangle into 
regions; we would like to label the regions by vertices of Gk , such that the labeling of the bottom part of T 
will correspond to q and the top part to q' , and then compute the matrix element T{T)qi^q from the labels. 
This is done as follows. 

The n marked points of a Kauffman n-diagram divide the top and bottom boundary into n+1 segments 
which we shall refer to as gaps. We shall say a set of gaps that bound the same region in the diagram are 
connected. For example, in the following Kauffman 3-diagram: 



7 6 5 4 




12 3 



the set of gaps {0,7} are connected as are the set of gaps {1,3,4,6}. We shall say that the pair {q',q) is 
compatible with T if once we label the gaps on the bottom from left to right by q{0), 5(1), • • • , q{n) and we 
label the gaps on the top from left to right by g'(0), q'{l), . . . , q'{n), then any set of connected gaps are all 
labeled by the same vertex of Gk. Thus in this case each region of T can be thought of as being labeled by 

a single vertex of Gfc. 

The matrix entry T{T)qi^g will only be nonzero in case the pair of paths {q',q) is compatible with T. 
In this case, the regions are indeed labeled by vertices in Gk; we can now do the following. To each local 
maximum and minimum of the Kauffman diagram T, we associate a complex number that depends on the 
labeling of the regions that surround them as follows: 



M l-l I I 

I I ' i+r ' i-r 

The matrix element T{T)qi,q at a compatible pair {q',q), is defined to be the product of the appropriate 
complex numbers over all local maxima and minima in T. 

For the map t(T) described above to be well-defined, it has to give the same result for isotopic Kauffman 
diagrams. An isotopic move can be seen to only create or eliminate local maxima and minima in pairs; we 
see that the conditions 

be+ice = 1 = ae-ide (1) 

are necessary and sufficient for the map to be isotopically invariant. A single extra constraint is needed to 
produce a representation of TL„(d) : 
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Claim 2.4 If the coefficients ae,be,ci,de satisfy and 

d = bfdi + ttiCi, (2) 
then the resulting map t defined as above is a representation ofTLn{d) . 

Proof: To prove the result we need only verify that the matrices T{Ei) satisfy the relations of Definition 
12.21 This amounts to verifying that the matrix elements of the operators on both sides of each relation are 
equal. Pictorially, a matrix element in a product of operators is given by stacking the operators together, and 
summing up over all possible labeling inside loops, such that the label inside the loop is different by exactly 
one from the label outside the loop. This summation corresponds to the summation over the intermediate 
index in matrix multiplication. 

We now check for the different relations in Definition 12. 21 For the generator Ei, there are only four types 
of non-zero elements, namely, four types of compatible pairs corresponding to the following types of labeling 











: l+l 


of the regions near the j-th strand: 


I 


I 


I 


I 




o ' 

l+l 


o ' 

l+l 







In the first two relations no loops are created when the operators are multiplied, so the verification follows 
from the isotopy invariance of r (i.e.. Equation (^). The third relation follows from Equation using the 
fact that one loop was created and there are two possible ways to label the region inside. □ 

We would like T{Ei) to be Hermitian, so that the induced representation on i?„ is unitary by Claim IT^ 
For this we add the constraints 

ae = Cg,bi = dg (3) 

Claim 2.5 // the coefficients ai, bi, c^, de satisfy Equations then T{Ei) — T{Ei)'^ for all i. 

Proof: We need to prove that T{Ei)q^qi = {T{Ei)q' _q)* . We need to check only for compatible pairs, namely 
for q, q' which on the i — 1, i, i + 1 entry are equal to ± 1, ^ for some i. This amounts to checking that 
a^di = b*f,c*g, bidi = b'^d^, and a^ci = a^c^. This follows from □ 
It is now left to solve Equations ^ - El to derive the definition of r. 

Claim 2.6 Define Xe = sin{TTe/k) for I e {1, fc - 1}. Then at ^ c} = ^ yTI' bi = d} = y^^^ satisfy 
Equations with d — 2cos{'K/k). 

Proof: We solve the equations directly. Define xi = j^^, and solve for xg. {xi is defined for £ G {1, . . . , fc — 
2}). We have -^^^ + xi — d, and if we set xq- --Xe — yi (and yi — whenever undefined) we have 
yi^2 + Ui — dyi-i- This is the familiar equation, defining the eigenvector of the adjacency matrix of the 
graph Gk 

fO 1 .. 
1 1 

1 

1 

\o ... 1 0/ 

with eigenvalue d. It is not difhcult to check that the (fc — l)-dimensional vector A is exactly this 
eigenvector. We have x^ = yt/yt-i,xi = a|, which gives — ^ X^j The Equations bg — l/ag_i and 
© determine bi,ci,di. □ 

Using these coefficients, we get the definition of T{Ei) as follows. T{Ei)q^q> = if {q, q') is not compatible 
with Ei. For a compatible pair {q,q'), T{Ei)q^q> is the product of two coefficients, one corresponding to the 
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maximum and the other to the minimum in Ei. For example, ii q = q' both move from the site £ to £ + 1 in 
the ith step, and then return to the site i in the i + 1th step, we have T{Ei)q^q = aicg. Given T{Ei), we can 
extend our definition of r to all elements in TLn{d). 

2.13 Unitary Path Model Representation of _B„ 

The previous section provided a representation r of TLn{d) defined for d — 2cos{'K/k). We have that 
d = —A^ — for A = ie~'"l'^^ , and also that r(i?j) are Hermitian. Thus, the conditions of Subsection 12. 71 
are satisfied. We define: 

Definition 2.13 The unitary path model representation of B„ is defined to be ^p{B) = t{pa{B)). 

The map can be extended to operate on tangles by letting pA be applied only to the crossings in the 
tangle. 

3 The Quantum Algorithm 

We are now ready to prove Theorems 11.11 and 11.21 Wc first translate the path model representation to work 
on qubits, and show that it can be implemented eSiciently. We use this to design the algorithms, and then 
prove their correctness. 

3.1 Moving to Qubits 

The adaptation of the path model representation to qubits is fairly straightforward. We simply switch from 
presenting paths by the list of their locations, to a binary representation which indicates the direction of 
each step. Thus, we shall interpret a string of n bits to be a sequence of instructions, where a shall mean 
take one step to the left and a 1 shall mean take one step to the right. We shall restrict our attention to 
those n bit strings that describe a path that starts at the leftmost vertex of Gk and remains inside Gu at 
each step. From here on when we say "path" we actually mean the bit string that represents the path. 

Definition 3.1 We define Pn,k,i to he the set of all paths p on Gk of n steps which start at the left most 
site and end at the £'s site. We define the subspace Ti.n,k,£ to be the span of \i) over all i e Pn.k.i- In a 
similar way, we define Pn,k to be all paths with no restriction on the final point, i.e., Pn.k — ^'i=iPn,k,i, and 
we define Tln,k to be the span of the corresponding computational basis states. 

We define a representation <1> as a homomorphism from TLn{d) to matrices operating on 'Hn,k- To define 
<I> it suffices to specify the images of the EiS, ^{Ei) ~ ^i. The operators $i are defined so that they 
correspond to the operators T{Ei) (see Subsection l2.12|) via the natural isomorphism between Vn,k and 'Hn,k- 
Thus, the transition from r to $ is merely a change of language. This is done as follows. 

To uniquely define $i on 7i„_fe, it suffices to define what it does to each basis element, namely, to \p) for 
P G Pn,k- We need the following notation: 

Definition 3.2 Let p\i denote the restriction of a path p to its first i — 1 coordinates. Given a path p on 
Gk, we denote by i{p) G {1, fc — 1} the location in Gk that the path p reached in its final site. Denote 
Zi = £ip\i). 

We can now define the operation of <I>i . <I>i is defined as an operation on the first i + I coordinates in a 
path p: 
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$ib|,oo) 

*i|p|zlO) 







(4) 



Azi+i 



A.. 
•\/ Azi+iA^ 

aT" 



-bl.oi) 







To apply <i>i on the rt-bit string \p) we tensor the above transformation with identity on the last n — i — \ 
qubits. For dealing with the edge cases, we use the convention Aj = for any j^{l,..,fc— 1}. 

Once we have defined $ is then extended to the entire algebra by the multiplication property of a 
representation, and by linearity. 

Claim 3.1 $ is a representation of TLn{d). are Hermitian. Moreover, $ induces a unitary representa- 
tion of the Braid group Bn, operating on 'Hn,k- 

Proof: The proof follows from the corresponding properties of t, namely Subsections 12.121 12.131 and the 
natural isomorphism between paths on Gk presented by their sequence of locations, Vn.fc, and paths presented 
as bit strings, namely Tin^k- ^ 

With the risk of confusion, we denote the unitary representation of i?„ induced by also by as in 
Definition l2.13l The only difference is that now ip{B) operates on TCn.k rather than on Vn,k- 



3.2 Efficient application of one crossing 

The matrices tpi are defined so far only on Tin.k which is a subspace of the Hilbert space of n qubits; we 
arbitrarily define their extension to the rest of the Hilbert space to be the identity. 

Claim 3.2 For all i e {l,...,n}, ipi can be implemented on the Hilbert space of n qubits using poly(n, k) 
gates. 

Proof: We note that the application of onp G Pn,k modifies only the i,i + l bits of p, and the modification 
depends on the location up to the ith step, namely on Zi — £{p\i-i). zi is a number which can be calculated 
efficiently and written on 0{log{k)) ancilla qubits, using the following standard technique: Initialize a counter 
register of, say, log{2k) qubits to the value 1. Then move along the qubits of the path p from left to right, 
and for each of the n qubits update the current state of the counter, i, by applying 

\b)\£) ^ \b)\£+{-lfmod2k) 

where b is the state of the currently read qubit. Since this is a unitary operation on log(2k) + 1 qubits, it can 
be applied using polynomially in k many elementary quantum gates (this is a standard result in quantum 
computation). We end up with the extra register carrying £{p\i^i). 

Now, ipi depends only on the location £{p\i-i) and on the i and i + I qubits. Hence, once again we have 
a unitary transformation which operates on logarithmically in k many qubits, and so we can implement it 
in polynomially in k many quantum gates. 

After we apply (pi, we erase the calculation of £{p\i-i) by applying the inverse of the first transformation 
which wrote the location down. □ 

As a corollary, we can deduce that 

Corollary 3.1 For every braid B e Bn, with m crossings, there exists a quantum circuit Q{B) that applies 
ip{B) on n qubits, using poly(m, n, k) elementary gates. 
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Proof: Order the crossings in the braid in topological order, and apply the corresponding unitary matrix of 
each crossing, tpi, one by one, in that order. Each crossing takes poly(n, /c) elementary gates by Claim rOl 
and there are m of them. □ 

3.3 The Algorithms 

We can now describe the algorithms. The input for both is a braid of n strands and m crossings, and an 
integer k. 

Algorithm Approximate-Jones- Plat-Closure 

• Repeat for j — 1 to poly{n, m, k): 

1. Generate the state la) = |1, 0, 1, 0, . . . , 1, 0) 

2. Output a random variable Xj whose expectation value is Tle{a\Q{B)\a) using the Hadamard test. 

• Do the same but for random variables j/j whose expectation value is Im{a\Q{B)\a) using the appro- 
priate variant of the Hadamard test. 

• Let r be the average over all Xj + iyj achieved this way. Output (— A)'^"'(-^'''^d^"/^^^Air/A^. 

Algorithm Approximate- Jones- Trace-Closure 

• Repeat for j = 1 to poly{n, m, k): 

1. Classically, pick a random path p e Pn^k with probability Pr(jp) cx A^, where £ is the index of the 
site which p ends at. 

2. Output a random variable Xj whose expectation value is TZ.e{p\Q{B)\p) using the Hadamard test. 

• Do the same but for random variables j/j whose expectation value is I'm{p\Q{B) \p) using the appropriate 
variant of the Hadamard test. 

• Let r be the average over all Xj + iyj . Output 

(_yl)3«.(B'-)^«-l^_ 



Claim 3.3 The above two quantum algorithms can be performed in time polynomial in n, m, k. 

Proof: Algorithm Approximate-Jones-Plat-Closure is clearly efficient, because the Hadamard test can be 
applied efficiently using Corollarv l3.1l To perform the first step of the second algorithm efficiently, we pick 
a random £ G {1, . . . , fc}, with probability proportional to A^, and then use Claim |^] below. We note that 
our calculations involve irrational numbers, Xi, but these can be approximated to within an exponentially 
good precision efhciently. □ 

Claim 3.4 Given n,k,£, there exists a classical probabilistic algorithm which runs in time polynomial in 
n and k, and outputs a random path in Pn,k,i according to a distribution which is exponentially close to 
uniform. 

Proof: (We thank O. Regev for a discussion that lead to this variant.) We briefly sketch the proof. We first 
note that the following nxk array 5*, defined by Sij — \Pi^k,j\ (the number of path on Gk of i steps that start 
at 1 and end at j) can be easily calculated efficiently using the recursion |-P„,fc,^| = l-Pn-i.fc.^-il + \Pn-i.k.i+i\- 
To pick a random path ending at £, we pick the values of the sites in p one by one in reverse order. If p{j) = £, 
we randomly decide if p{j — 1)=^— lori!-|-l according to the ratio S^-i^^-i : S'j-i^^+i. Given p{j — 1) we 
can continue to pick p{j — 2) in a similar manner, and so on. □ 
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Theorem 3.1 ^Z^orii/im Approximate-Jones-Trace-Closure approximates the Jones polynomial of B*'^ at A ^ = 
^2-iTi/k ^ within the precision specified in Theorem ] 1. 11 

Proof: We will need the following definition. 

Definition 3.3 Define Tr„(VF) for every W in the image of ^{TLn{d)) to be: 

^ fe-i 

Tr^{W)^-Y,^iTr{WU) 

e=i 

where W\i denotes the restriction ofW to the subspace 'Hn,kj., md Tr denotes the standard trace on matrices. 
The renormalization is N — '^iXidim{TLn,k,i) where the sum is taken over all (.'s such that Pn,k,e is non 
empty. 

This definition makes sense because matrices in the image of Lp are block diagonal, with the blocks indexed 
by the last site of the paths: 

Claim 3.5 For any T e TL^id), <i>{T)nrM C W„,fc,^. 

Proof: <i>i cannot change the final point of a path since it only moves 01 to 10 and vice versa. □ 

Hence, the above trace function simply gives different weights to these blocks (and gives zero weights on 
strings that aren't paths). We claim that Tr„ is a Markov trace. 

Claim 3.6 The function Trn(-) satisfies the three properties in Claim [^31 

Proof: That rr„(<l>(l)) = 1 follows from the renormalization. The second property follows from Claim 1^31 
plus the fact that the standard trace on matrices satisfies this property, so Tr„(-) satisfies it on each block 
separately. To show the Markov property, we have to show that if X e TLn-i{d)) then rr„((f>(X)(f>(ii'„_i)) = 
^Tr„($(X)). We note that for any X G Ti„_i((i), ^{X) can be written as a linear combination of terms 
of the form \p){p'\ (8i /, with p,p' £ Pn~i,k and the identity operates on the last qubit. By hnearity, it 
suffices to prove the Markov property on such matrices. Writing \p){p'\ ^ I = \pO){p'0\ + \pl){p'l\ we 
require: Tr„(|pO) (p'0|$„_i + |pl) (p'l|$„_i) = ^rr„(|pO) (p'0| + \pl)(p'l\). This can be easily verified using 
the definition of $ by checking the two cases p — p' and p ^ p' . 

We start with the case p ^ p'. In this case the right hand side is 0. As for the left hand side, {p'0\^n-i 
has a zero component on {pO\. To see this, we check the two cases: if p' ends with then {p'0\^n-i — 0, 
otherwise p' ends with 1. (p'0|$„_i is then a sum of two terms, one equals to {p'0\ and is therefore different 
than (pO|, and the other ends with 01 and is thus also different from {pO\. The same argument works to show 
that (p'l|$„_i has a zero component on {pl\. Hence, the left hand side is also 0. 

It is left to check the equality in the case p = p' . We require 

rr„(|pO)(pO|<f„_i + bl>(pl|<i>„_i) - lTr„(|pO)(pO| + bO>(pO|). 

d 

Suppose £{p) = £. Then the right hand side is equal to ^(Af_i + A^-i-i) — A^, using the properties of the 
eigenvector A, as in Claim [TEl To see that the left hand side is the same, we again divide to cases. Suppose 
first that p ends with 0. In this case (pO|$„_i = 0. As for the other term, (pl|$„_i — j^(pl|, using the 
definition of ^n-i and the fact that p without its last step ends in ^ -f 1. The weight in the trace of the left 
hand side is A^+i, and so the left hand side is equal to A^ too. The argument is similar in the case that p 
ends with 1. □ 

By the uniqueness of the Markov trace. Lemma [2. II we have that Trn{(p{B)) — tr{pA{B)). Hence, using 
Lemma [2.21 we have that for any braid B £ B^ 

Lemma 3.1 
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Due to Lemma T-j. II the correctness of the algorithm follows trivially from the following claim. 

Claim 3.7 With all but exponentially small probability, the output r satisfies \r — Tr„((/3(i?))| < e for e 
which is inverse polynomial in n, k,m. 

Proof: The Hadamard test indeed implies that the expectation of Xj for a fixed p is exactly Tle{p\Lp{B)\p) . 
The expectation of the variable Xj taken over a random p is thus 

Ei,peP^,,^,>^enei{p\cp{B)\p)) _ J2,Xine{Tri^iB)U)) _ 
= ne{Tr,MB)))- 

The same argument works for the imaginary part. Since r is the sum of two averages of polynomially many 
i.i.d random variables, each taking values between 1 and —1, the result follows by the Chernoff-Hoeffding 
bound. □ 

This completes the proof of Theorem 13. II □ 

Theorem 3.2 ^/^orif/im Approximate-Jones-Plat-Closure approximates the Jones polynomial of B^^ at A^^ = 

g27ri/fe^ within the desired precision as in Theorem M.^ 

Proof: By the correctness of the Hadamard test, and by the Chernoff-Hoeffding bound, the variable r which 
the algorithm computes is, with exponentially good confidence, within 6 — l/poly{n, m, k) from {a\ip{B)\a) , 
which is equal to Tr{Lp{B)\a) {a\) . We need to connect this expression to the Jones polynomial of the plat 
closure of B. The main observation here is that the plat closure of a braid B is isotopic to the trace closure 
of a tangle C achieved by applying the braid on n/2 capcups, as in the following picture: 




It therefore suffices to relate the Jones polynomial of C*'' to Tr{(f{B)\a){a\). Since the question is now 
stated in terms of trace closures and traces, there is hope to be able to apply Lemma f2. 21 and Claim IT^ as 
in the proof of Theorem 13. II But we first need to make the connection between the projection on \a) and 
capcups. 

Claim 3.8 \a){a\ = $i$3 . . .$„_i/d"/2. 

Proof: It is easy to verify that $i$3 . . . $„-i applied to any path except for \a) gives 0, and when applied 
to I a) it gives the desired factor. To do this we use the fact that $i commute if their indices are more than 
one apart, and so wc can first apply <I>i, then <I>3 and so on. Since the path starts at the left most site, $i 
on p simply applies the following rescaled projection: d|10)(10| on the first two coordinates. This projection 
forces the first two coordinates to be 10, and so ^i\p) returns to the starting point after two steps. Therefore 
a similar argument applies when we apply $3 on the next two coordinates, and so on. By induction, we get 
the desired result. □ 

We thus have, using Definition 13.31 and Claim 

N 

{a\^{B)\a) ^Tr{^{B)\a){a\) - -Tr^{^{B)\a) {a\) = 
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N 



rr„((p(B)$i$3...$n-i/d"/') 



N 



Aid"/2 



By the uniqueness of the Markov trace, Lemma l2.ll and by Claim 1X1)1 we have that Tr„(<p(C)) — 
tr{pA{C)). Using Lemma [2. 21 we have: 



which we substitute in the previous equation. This completes the proof of Theorem 13.21 □ 
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